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Search for double charmonium decays of the P-wave spin-triplet bottomonium states 
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The order of magnitude discrepancy between the cross 
sections of the double charmonium production processes 
e+e" -^ J/ipr]c, J/^ic±Jp^c, i^'Vc, J/i^XeO, and ■(/''Xco 
measured in the Belle |l|, |2| and BaBar 3] experiments 
and those of the leading-order non-relativistic QCD 
(NRQCD) predictions [^-Ig has been a great challenge 
to theorists. After great efforts it was shown that agree- 
ment can be achieved by taking into account QCD radia- 
tive and relativistic corrections [J, |7H1C| . 

As in e~^e~ annihilation, double charmonium final 
states can also be formed in bottomonium decays, which 
provide a new test of the dynamics of hard exclusive 
processes and charmonium structure. While the rate 
of rji, — ^ J/ipJ/ip decay was calculated a long time ago 
by many authors 11], the rates for P-wave spin-triplet 
bottomonium states Xbj {J=0, 1, 2) decays into double 
charmonium states were calculated only recently using 
various theoretical models after a first attem pt a bout 30 
years ago with a perturbative QCD method [l2[. 

The authors of Ref. [13'| calculated XbJ -^ J/i^J/i' in 
the framework of the NRQCD factorization formalism, 
including second-order relativistic corrections in the rel- 
ative charm-quark velocity Vc, as well as an electromag- 



netic correction. The branching fraction is predicted to 
be of order 10"^ for Xfco or x&2 — > J ji^J 1"^, and 10~^^ 
for Xh\ —> J/ipJ/ip. The authors of Ref. [14] considered 
corrections to all orders in Vc in the charmonium rest 
frame, and found decay partial widths that are about 
a factor of three larger than those in Ref. 13j. In the 
light cone (LC) formalism, however, much larger produc- 
tion rates (with uncertainties of 50% to 70%) are ob- 
tained in Ref. ^■. B(xb.j ^ J/^J/^) = 9.6 x IQ-^ or 
1.1 X 10-3, B{xbj -^ Jl^^') = 1-6 X lO-'^ or 1.6 x lO^^, 
and B{xhJ -^ V-V) = 6.6 x IQ-^ or 5.9 x lO"'* for J=0 
or 2, respectively. These results are not very different 
from the perturbative QCD calculation with the relative 
motion of the charm-quark in the Xb decays taken into 
account [16|. It is therefore necessary to pin down the 
source of such a significant difference between various 
models. 

In this paper, we report a search for Xbj decays to 
double charmonium in T{2S) radiative transitions, i.e., 
T{2S) -^ -fXbj -> jJ/ipJ/ip, iJ/^i'', and 7V'V. In 
order to detect the signal efficiently, we only require 
one J/ijj or -0' candidate to be fully reconstructed (or 
"tagged"). We require that the missing mass of the 



J /ip (or ■)/)') and the radiative photon candidate be in 
the J /^ or if)' mass region. The missing mass is de- 
fined as Mmiss = \J{Pe+e- — PfY^ where Pe+e- is the 
4-momentum of the e+e~ colhsion system, and Pf is the 
sum of the 4-momentum of the observed final-state par- 
ticles. Double counting of an event is allowed if both 
charmonium states satisfy the tag criteria. The proba- 
bility of double counting depends on the final states and 
varies from a few per mille to less than 5%; this is taken 
into account in the efficiency estimation using MC sam- 
ples generated with both charmonium states decaying 
generically. For jj/ipj/ip, we reconstruct one J/ip signal 
from £^£~ (£ = e or /z) and require the missing mass of 
'jJ/ip be within the J/ip mass region. For ^J/ipip', three 
modes are included: (1) J/V" — >■ ^"'"€~ with Mmiss(7J/V') 
required to be within the ip' mass region; (2) xp' — > 
TT+TT" J/V' -^ n+TT-i+i- with MmiBsil^') withlu the J/ip 
mass region; (3) ip' — ?> £^i^ with M^-nissi'yip') within 
the J/ip mass region. For ^ip'ip', two modes are used: 

(1) one ip' — >■ TT^TT~J/ip -^ TT^TT^£^£^ is identified and 
-^miss(7V'') is required to be within the ip' mass region; 

(2) one Ip' — >■ £'^£~ is identified and M^issiji/j') is re- 
quired to be within the ip' mass region. 

This analysis is based on a 24.7 fb'^ T(2S') data sam- 
ple (158 million T(2S') events ^), and a 89.4 fb^i con- 
tinuum data sample collected at y/s = 10.52 GeV. Here 
y/s is the center-of-mass (CM.) energy of the colliding 
e+e". The data are collected with the Belle detector [18| 
operating at the KEKB asymmetric-energy e+e" col- 
lider [13]. 

EVTGEN [20| was used to generate Monte Carlo (MC) 
simulation events. For signal MC samples, the angular 
distribution for T{2S) — > "fXbJ is simulated assuming 
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a pure El transition {dN / d cos 9 j 

1, — i, ^ for J = 0, 1, 2, respectively [21 
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Here 9^ 
is the polar angle of the T(25') radiative photon in the 
e"'"e~ CM. frame. Uniform phase space is used for Xbj 
decays [22 1. 

Since no experimental measurements are available for 
the widths of Xbj |2J], and the theoretical expectations 
are at 1 MeV level or less 13| , the widths of XbJ are set 
to be zero. Generic decay modes are used for the J/ip 
and Ip'. T(2S) MC events with generic decays produced 
with PYTHIA [25| with two times the effective luminosity 
of data are used to check the possible backgrounds from 
T(2S') decays. 

The detector is described in detail elsewhere [l8|. It 
is a large-solid-angle magnetic spectrometer that consists 
of a silicon vertex detector (SVD), a 50- layer central drift 
chamber (CDC), an array of aerogel threshold Cherenkov 
counters (ACC), a barrel- like arrangement of time-of- 
fiight scintillation counters (TOF), and an electromag- 
netic calorimeter comprised of CsI(Tl) crystals (ECL) 
located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux-return located 
outside the coil is instrumented to detect K^ mesons and 



to identify muons (KLM). 

For well reconstructed charged tracks, the impact pa- 
rameters perpendicular to and along the beam direction 
with respect to the nominal interaction point are required 
to be less than 0.5 cm and 4 cm, respectively, and the 
transverse momentum in the laboratory frame is required 
to be larger than 0.1 GeV/c. We require the number 
of well reconstructed charged tracks to be greater than 
three for ^ J/ip J/ip, and greater than four for ^J/ipip' and 
^ip'ip' . For the modes with ip' in the final states, events 
with exactly four charged tracks are removed to suppress 
the significant background from QED processes. For each 
charged track, information from different detector sub- 
systems is combined to form a likelihood Ci for each par- 
ticle species [26|. A track with TZk = 



Ck+C^ 



< 0.4 is 



identified as a pion with an efficiency of about 97% for the 
momentum range of interest; about 3.5% are misidenti- 
fied K tracks. For electron identification, the likelihood 



ratio is defined as TZe 



c^+c^ 



where £e and L^ are the 



likelihoods for electron and non-electron, respectively, de- 
termined using the ratio of the energy deposit in the ECL 
to the momentum measured in the SVD and CDC, the 
shower shape in the ECL, matching between the position 
of the charged track trajectory and the cluster position 
in the ECL, the hit information from the ACC and the 
dE/dx information in the CDC [27]. For muon identifica- 
tion, the likelihood ratio is defined as 7?.^ — _^ _|_£^i£ , 
where >C^, C-j^, and Ck are the likelihoods for muon, pion, 
and kaon hypotheses, respectively, based on the match- 
ing quality and penetration depth of associated hits in 
the KLM [23. 

A neutral cluster is used as a photon candidate if it 
docs not match the extrapolation of any charged track 
and its energy is greater than 50 MeV. In calculating the 
recoil mass of ^J/ip or ^ip' , all photon candidates except 
those within 0.05 radians of the electron/positron tracks 
are included. No tt'^ signal is observed in combining the 
low energy radiative photon with any of the remaining 
photon candidates in the event after all the selection cri- 
teria are applied. 

In order to correct for the effect of bremsstrahlung 
and final-state radiation, photons detected in the ECL 
within 0.05 radians of the original e"*" or e~ direction 
are included in the calculation of the e+/e^ momentum. 
For the lepton pair used to reconstruct J /ip, both tracks 
should have TZe > 0.95 in the e+e" mode; or one track 
should have Tl^ > 0.95 while the other should satisfy 
TZ^ > 0.05 in the ^l^ ^i~ mode. The lepton pair identifi- 
cation efficiency is about 90% for J/ip — >■ e+e" and 87% 
for J/i]j — )• //+/i^. In order to improve the J/ip momen- 
tum resolution, a mass-constrained fit is then performed 
for J/ip signals in all the modes. As different modes 
have almost the same J/ip mass resolutions, the J/ip sig- 



nal region is defined as \MgJrg~ 



T-j/ii 



< 0.03 GeV/ 



(w 2.5cr), where mj/^ is the nominal mass of J/ip 



The J/t/i mass sidebands are defined as 2.97 GeV/c^ < 
Mt+E- < 3.03 GeV/c^ or 3.17 GeV/c^ < M^+£- < 
3.23 GeV/c^, and are twice as wide as the signal re- 
gion. For ip' — )■ €+£~, the -0' signal region is defined 
as \Mi+p- - m,0/| < 0.0375 GeV/c^ (« 2.5cr), where m^- 
is the nominal mass of V'' |2^- The ip' mass sidebands 
are defined as 3.535 GeV/c^ < M^+^- < 3.610 GeV/c^ or 
3.760 GeV/c2 < M^+i- < 3.835 GeV/c^, and are twice 
as wide as the signal region. For ij^' — > tt+tt^ J/'0, we 
require the two pion candidates be positively identified. 
The if)' signal region is defined as iM^+jr-j/^ — m^'| < 
0.009 GeV/c^ (« 3cr). Figure [1] shows the mass dis- 
tributions of the reconstructed J/^ — )■ £+£~ (a), V' ~^ 
7r+7r~ J/?/) (b) and ^/i' — >■ £+£~ (c) candidates. 

Figure [2] shows scatter plots of the photon spectra 
in the e+e^ C.M. frame versus (a) Mmiss{lJ/i') with 
J/ip -> ^+£" reconstructed, (b) M^issil"^') with -0' -^ 
TT+TT" J/'0 reconstructed, and (c) Mmiss(7V'') with tp' — >• 
£+£~ reconstructed. No evidence for J/ip or ip' signals 
can be seen in the ^J/ip or 7?/;' missing mass distribu- 
tions. 

Figure |3] shows the simulated photon spectra in the 
e+e~ C.M. frame from the T(25') -> 7Xbj -> ^JjipJlip 
MC samples. Breit-Wigner (BW) functions convolved 
with Novosibirsk functions |29| are used as XbJ signal 
shapes while Chebychev polynomial functions model the 
combinatorial backgrounds (^ 11% in Xhj signal re- 
gion). The extended maximum likelihood fits to the 
photon spectra with all the parameters free are shown 
in Fig. [H Based on the fit results, the efficiencies are 
(5.75 ± 0.12)%, (6.25 ± 0.12)%, and (5.87 ± 0.12)% for 
T(2S') -> 7Xbj -^ ^JjipJlip for J = 0, 1 and 2, respec- 
tively. Similarly, the sum of the efficiencies from all the 
modes is found to be (3.40 ± 0.06)%, (3.78 ± 0.06)%, 
and (3.53 ± 0.06)% for T(2S') ^ 7X6j ^ iJ/ip-ip', 
(2.06 ± 0.04)%, (2.15 ± 0.04)%, and (2.09 ± 0.04)% for 
T(2S') — > 7XhJ — > 7V'V for J = 0, 1 and 2, respectively. 

After all the event selections, no events from the T(25') 
MC sample with generic decays survive. Other possible 
backgrounds with J/0/-0' signals from channels such as 
e+e" — >■ J/ipXcJ, ip'Kr.u have very small cross-sections 
(at the few fb level [30[) and hence are neglected in the 
analysis. 

Figures EKa), (b) and (c) show the photon spectra from 
T(2S') data for Xbj -^ J/tpJ/tp, J/ipip', and ip'ip' candi- 
date events, respectively, with all the modes included. 
Here the shaded histograms show the J/-0 or -0' mass 
sidebands normalized to the width of the J/ip or ip' signal 
range and the dashed histograms are the normalized con- 
tinuum contributions. The continuum background con- 
tribution is extrapolated down to the T(2S') resonance. 
For the extrapolation, three factors are applied to ac- 
count for: the relative luminosities of the two samples, 
efficiency dependence on the C.M. energy, and cross sec- 
tion dependence on the C.M. energy. The cross section 
extrapolation with C.M. energy is assumed to have a 1/s 



dependence. 

No clear XbJ signals are observed in Fig. U) For 
XbJ -^ Jl'^Jl'4', a unbinned extended maximum like- 
lihood method is applied to the photon spectrum with 
the MC simulated signal shape smeared with a Gaus- 
sian function to take into account a 8.5% difference in 
photon energy resolution between data and MC sam- 
ples. The photon energy resolution is measured with 
T(2S') -^ -iXbj -)■ 77T(1S'), T(15) -^ ii+iJ.- events. For 
XbJ -^ J /ipip' {ip''4'') decays, an unbinned extended max- 
imum likelihood simultaneous fit is performed to all the 
modes mentioned above. The ratios of the XbJ yields in 
different modes are fixed to e^ (i denotes the i-th mode) 
with all the intermediate state branching fractions in- 
cluded, and £i is the MC-determined efficiency for the 
i-th mode. The fits are performed with the same method 
as in the J /ipJ /"^P mode. Figure U shows the fit results, 
where for (b) and (c) the solid curves are the sum of 
all the fit contributions, and the dashed curves are the 
sum of the background functions. In all of the modes, the 
background levels from the fits are a little higher than the 
estimations from the normalized continuum or the nor- 
malized J /ip /^' mass sidebands. It may indicate that 
there are double-charmonium production together with 
one photon or more particles in T(2S') decays. 

The upper limit on the number of signal events at 
the 90% C.L. (n"P) is calculated by solving the equa- 

/ C{x)dx 

tion °_foo =0.9, where x is the number of signal 

J„ C(x)dx 

events, and C{x) is the likelihood function depending on 
X from the fit to the data, with x being the number of 
signal events in the fit. The values of n"P are found to 
be 21, 13, and 22 for Xbj ^ J/^J/V'; 20, 5.8, and 17 for 
XbJ — ^ J/i'tp'', and 3.0, 12, and 3.3 for XbJ -^ V'Vj for 
J — 0, 1, and 2, respectively, when requiring the signal 
yields to be non-negative in the fit. 

There are several sources of systematic errors for the 
branching fraction measurement. The uncertainty in the 
tracking efficiency for tracks with angles and momenta 
characteristic of signal events is about 0.35% per track, 
and is additive. The photon reconstruction contributes 
an additional 3.8% per photon. The uncertainty due to 
particle identification efficiency is 1.3% for each pion in 
-0' -^ TT^TT^J/ip. According to a measurement of the 
lepton identification efficiency using a control sample of 
77 — >■ £'^£~, the MC simulates data within 1.7% for an 
electron-positron pair and 1.7% for a muon pair. Ac- 
cording to MC simulation, the trigger efficiency is greater 
than 99.5% and we take 0.5% as systematic error due to 
the trigger simulation uncertainty. Errors on the branch- 
ing fractions of the intermediate states are taken from 
the PDG "gj], which are about 12%, 6.0% and 5.0% for 
Xbo, Xbi and Xb2 decays. By changing the order of the 
background polynomial and the range of the fit, the rela- 
tive difference in the upper limits of the number of signal 
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FIG. 1: The reconstructed J/ip — >■ ^"""^ (a), ip' ~^ tt'^tt J/ip (b), and tp' — > ^"'"Z (c) candidates mass distributions from data. 
The arrows show the required signal mass regions. 
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FIG. 2: Scatter plots of the photon spectra in the e^e~ C.M. frame versus (a) Mmiaai'yJ/tp) with J/tp — >■ £'^£~ reconstructed, 
(b) Mmiss(7^') with ip' — > -k^tt' J/ip reconstructed, and (c) Mmiss(7^') with tp' — > i^£~ reconstructed. The dotted lines show 
the J/^p or ?/)' signal regions (~ ±3(t). 



events is 7.6%-28% depending on the decay mode, which 
is taken as systematic error due to the uncertainty of fit. 
For our MC signal samples, J/4' and V'' decays are simu- 
lated with a generic decay model. The signal efhciencies 
are determined based on the fitted results. The error 
on the number of fitted signal events is less than 2.1%, 
which is taken as the MC statistical error in the efficiency. 
The masses of XbJ have been measured well 24| and the 
uncertainties on the masses of Xbj do not affect the effi- 
ciency determination. Comparing several theoretical cal- 
culations, the maximum values of Xfco and Xb2 widths 
are 2.15 MeV/c^ [31] and 0.33 MeV/c^ [32l], respectively. 
The efficiency differences between these values and the 
nominal values are taken as systematic errors due to the 
uncertainty of resonance parameters, which are less than 
5.2% and 1.2% for Xco and Xb2 decays. Finally, the un- 
certainty on the total number of T(2S') events is 2.3%. 
Assuming that all of these systematic error sources are 
independent, and combining them in quadrature, we ob- 
tain the total systematic error listed in Table HI 

Since there is no evidence for signals in the modes stud- 
ied, we determine upper limits on the branching fractions 
of XbJ to double charmonia. Table H] lists the upper lim- 
its n"P for the numbers of the signal events, detection 
efficiencies, systematic errors, and upper limits on the 
branching fractions of Xbj decays. In order to calculate 
conservative upper limits on these branching fractions, 
the efficiencies are lowered by a factor of 1 — cTgys in the 
calculation. 

To summarize, we find no significant signals in the 



XbJ ~^ J/i^J/i^i J/i^ip' : or V'V final states using a sam- 
ple of 158 million T(2S') events. The results obtained on 
the XbJ decay branching fractions are listed in Table ID 
Our upper limits are much lower than the central values 
predicted in the LC formalism 15] and pQCD calcula- 
tion [l6|, but are consistent with calculations using the 
NRQCD factorization approach Ij, [l4 1 . 



TABLE I: Summary of the limits on Xbj decays into J/%pJ/%p, 
J/^ptp' , and ^p'tp' . Here n"'' is the upper limit on the number 
of signal events, e is the sum of the efficiencies from different 
modes with J/^p and ^p' decay branching fractions and trigger 
efficiency included, asys is the total systematic error, and Br is 
the upper fimit on the branching fraction of Xbj decays, where 
the values of B{T{2S) -^ jxbj) = (3.8 ± 0.4)%, (6.9 ± 0.4)% 
and (7.15 ± 0.35)% for J = 0, 1 and 2 are used ;24]. The 
upper limits are at 90% C.L. 



Channel 


n"P 


e(%) 


asys(%) 


Br 


Xm) -> J/ipJ/ip 


21 


5.8 


16 


7.1 X 10"'' 


Xbi -> J/ipJ/ip 


13 


6.3 


30 


2.7 X lO-'^ 


Xb2 -^ J/ipJ/i^ 


22 


5.9 


27 


4.5 X lO"'' 


XM) -^ J/Tp1p' 


20 


3.4 


17 


1.2 X 10"^ 


Xbl -^ J/i>1p' 


5.8 


3.8 


15 


1.7 X 10"^ 


Xb2 -S- J/lp1p' 


17 


3.5 


16 


4.9 X lO-'^ 


XbO -i> l/'V 


3.0 


2.1 


20 


3.1 X 10"^ 


Xbl -^ ip'ip' 


12 


2.2 


17 


6.2 X lO-'^ 


X62 -i> 'ip'tp' 


3.3 


2.1 


12 


1.6 X 10"^ 



We thank the KEKB group for excellent operation 
of the accelerator, the KEK cryogenics group for effi- 










> 1000 


- (^) ft 


MC 


o 


1 1 




S 800 


f 




CM 


j 




« 600 


f 










1 400 


- / ' 




> 


f I 




"" 200 

n 


^ L, 





0.08 0.1 0.120.140.160.18 0.2 
Ev(GeV) 



0.08 0.1 0.120.140.160.18 0.2 
E, (GeV) 



0.08 0.1 0.120.140.160.18 0.2 
Ey(GeV) 



FIG. 3: The fits to the photon spectra from the T{2S) — > "fxtj — >■ 'yJ/i^J/ip MC signal samples with one J/tp reconstructed 
and the yJ/tp recoil mass within the J/ip mass region for (a) Xbo, (b) Xbi and (c) Xb2, respectively. The XbJ shapes are 
described by Breit-Wigners convolved with Novosibirsk functions, while Chebychev polynomial functions are used to describe 
the background. 
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and dashed histograms are normalized continuum contribu- 
tions. The fits to the photon spectra are described in the 
text. The solid curves are the best fits, the dashed curves 
represent the backgrounds. The arrows show the expected 
central positions of the Xbj states. 
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